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ABSTRACT 


This  report  presents  the  results  of  studies  on  the  stiffness  variation  effects  to  the 
structures  such  as  ship,  missile  and  airplane  subjected  to  an  axial  or  a  vertical  impact  load. 
The  theoretical  studies  were  carried  out  for  a  simply  supported  beam  under  an  axial 
compressive  load  and  a  cantilever  beam  under  a  vertical  load.  A  transient  analyses  for 
beams  and  cylindrical  shells  were  performed,  and  the  results  were  evaluated.  It  was  also 
investigated  the  effect  of  frequency  variation  of  an  impact  load  on  the  response  of 
structures.  A  finite  element  modeling  and  a  numerical  analysis  were  performed  by 
MSC/PATRAN  and  MSC/NASTRAN  program. 
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I.  INTRODUCTION 


Dynamic  buckling  of  structures  has  been  studied  in  the  field  of  naval  architectural, 
aerospace,  nuclear  power  plant  and  civil  engineering.  Transient  or  impulsive  loads,  which 
may  exceed  the  critical  static  Euler  limit,  can  also  be  an  important  design  consideration. 
Especially,  in  case  of  warship  or  submarine,  the  design  of  structures  to  resist  shock  loads 
from  hostile  weapons  or  missile  launching  is  important  for  the  protection  of  crew's  safety 
and  environment.  On  the  ground  of  that,  many  people  have  researched  dynamic  buckling  of 
structures. 

Takuo  and  Yukio  [Ref.  1]  studied  a  simply  supported  column  with  an  initial 
deflection,  which  had  one  end  struck  by  a  mass  along  the  axial  line  at  higher  impact 
velocities. 

Tanchum  and  Haim  [Ref.  2]  calculated  and  compared  the  buckling  problems  of 
beam  and  plate  subjected  to  an  axial  impact.  They  studied  on  the  initial  geometric 
imperfection  of  structures  and  determined  a  dynamic  load  amplification  factor  by  analytical 
studies  with  ADINA  program. 

N.  G.  Pegg  [Ref.  3]  investigated  the  effect  of  a  ring  stiffener  and  its  size  and  spacing 
on  the  dynamic  buckling  response  of  a  cylinder. 

W.  Gu  and  W.  Tang  [Ref.  4]  studied  the  dynamic  plastic  buckling  of  cylindrical 
shells  subjected  to  general  external  impulsive  velocity  and  subjected  to  asymmetric 
impulsive  loadings. 

Murli  and  Judah  [Ref.  5]  investigated  the  dynamic  buckling  of  geometrically 
imperfect  columns  with  viscous  damping  under  an  axial  compressive  pulse.  Viscous 
damping  effects  became  significant  for  extremely  short  pulses  and  that  higher  damping 
caused  an  increase  in  the  dynamic  buckling  load  and  a  smoother  deflection  pattern. 


1 


S.  Kenny  and  N.  Pegg  [Ref.  6]  did  the  paper  concerned  with  dynamic  elastic  pulse 
buckling  events,  characterized  by  the  unacceptable  growth  of  lateral  displacements  due  to 
an  intense  transient  load  that  significantly  exceeds  the  critical  static  Euler  limit. 

O.  Aksogan  and  A.H.  Sofiyev  [Ref.  7]  studied  the  dynamic  buckling  of  an  elastic 
cylindrical  shell  with  variable  thickness  subjected  to  a  unifonn  external  pressure  which  was 
a  power  function  of  time. 

Shijie  and  Hong  [Ref.  8]  investigated  the  dynamic  elastic  buckling  of  simply 
supported  columns  subjected  to  intermediate  velocity  impact.  They  proposed  the  dynamic 
buckling  criterion  to  determine  the  critical  buckling  condition  and  to  estimate  the  dynamic 
buckling  critical  load. 

Galib  H.  and  Christos  C.  [Ref.  9]  suggested  the  computational  simulation  method  to 
evaluate  the  detenninistic  and  non-deterministic  dynamic  buckling  of  adaptive  composite 
shells. 

As  mentioned  above,  there  are  many  studies  for  a  continuous  beam  or  a  cylinder 
under  axial  shock  load.  But  it  is  not  easy  to  find  the  paper  on  the  stiffness  variation,  one  of 
important  design  factors  in  the  preliminary  design  stage  of  the  structures  such  as  ship, 
airplane  and  missile.  Generally,  the  body  section  stiffness  of  those  are  different  each  other 
by  access  doors,  stiffener,  skin  cutout  and  so  on.  Therefore,  the  cost  by  a  trial  and  error  can 
be  reduced  if  the  structure  design  process  considers  the  stiffness  distribution  of  a  system  in 
preliminary  design  step.  In  this  report,  we  investigated  the  effect  of  the  stiffness  variation 
on  the  stepped  slender  column  subjected  to  an  axial  or  a  vertical  impact  load.  The 
theoretical  studies,  for  a  simply  supported  beam  under  an  axial  load  and  a  cantilever  beam 
under  a  vertical  load,  were  carried  out.  A  transient  analyses  for  beams  and  cylindrical  shells 
were  done,  and  the  results  were  evaluated.  Also,  it  was  investigated  the  effect  of  the 
frequency  variation  of  an  impact  load  on  the  response  of  structures.  A  finite  element 
modeling  and  a  numerical  analysis  were  performed  by  MSC/PATRAN  and 
MSC/NASTRAN  program. 
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II.  MATHEMATICAL  FORMULATION 


This  chapter  is  concerned  with  the  equations  of  motion  of  a  stepped  beam  subjected 
to  axial  or  vertical  loads.  Mathematical  formulations  are  carried  out  for  a  simply  supported 
beam  and  a  cantilever  beam  as  follows. 


A.  A  SIMPLY  SUPPORTED  BEAM 

From  H.E.  Lindberg  and  A.L.  Florence  [Ref.  10],  the  equations  of  motion  of  a 
simply  supported  beam  under  an  axial  compressive  force  P  are 

d4v  d2  d2v 

£/a?+F  a?('v+^)  +  ^aF^°  (1) 

Where,  E  is  an  elastic  coefficient,  I  is  an  area  moment  of  inertia  and  v0  is  an  initial 
deflection,  p  is  a  mass  density  and  A  is  an  area  of  beam. 

For  static  buckling,  the  inertia  term  is  neglected  and  Equation  (1)  becomes 

eiAz+pAz=.pAisl 

dx  dx  dx~ 

,  P 

Or,  substituting  k  =  — 

d'y  |  k 2  d2y  _  k2  d2y0 
dx 4  dx2  dx2 

Homogeneous  solution  is 

yh  =  Cj  cos  kx  +  C2  sin  kx  +  C3x  +  C4  (4) 

And,  we  can  get  Cj ,  C2 ,  C3  and  C4  from  the  boundary  conditions. 

In  case  of  the  stepped  beam  as  in  Figure  1,  a  homogeneous  solution  is 

3 

yh  =  cos  kjX  +  Cj2  sin  kjX  +  Cj3x  +  C jA  ]  (5) 

j= i 


(2) 


(3) 
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Figure  1.  Stepped  Simply  Supported  Beam  Subjected  to  Compressive  Force 


Applying  the  boundary  and  the  continuity  conditions  to  get  the  coefficients  of  Equation  (5) 
At  x  =  0 


v,(0)=0  (6) 

=  0  (7) 

OX 

At  x  =  xt,x2 

y,(x1)  =  y2(xl)  (8) 

y[(xl)  =  y'2(xl)  (9) 

y''(x])  =  y"(x])  (10) 

y;"(x1)  =  y'"(x1)  (11) 

v2(x2)=y3(x2)  (12) 

y2(x2)  =  yj(x2)  (13) 

yi(x2) =y"(x2)  (14) 

y2(x2)=y"'(x2)  (15) 

At  x  =  x2 

y3(x3)  =  o  (16) 

It  y,M  =  o  07) 
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Rearranging  Equations  (5  through  17)  by  the  matrix  fonn  as  follows, 
h,M2]iC*M°}  j=l-2,3  and k  =  1,2, 3, 4 

And  are 

[?ilOl]=l.  [^0104 ]=  1  ’  [^020 1  ] =  K  ’ 

[-^0302  ]=  sin  ktXt  ,  [^0303  ]=  'fi  ’  [^0304  ]=  ^  ’ 


[^0407  ] 


[^0708  ]  =  1  ’ 

[^0712  ]=  1  ’ 

[7o809  ] — k2  sin  ky 


r  ]--k2  \: 

0201 J  "-1  5  LJ 

^0301  ]=cosk1x1 , 

[■^0303  ]  =  X\  ’  D 

r  l-i  f 

0304  J  1  5  V 

^0305  ]=  ~~  COS  k2xx , 

[^0307  ]=  ~~Xt  ’  [J 

^0308  ]~  1  ’  [ 

^0401  ]  ^1  sinkj.Tj , 

[^0403  ]=^’ 

[^0405  ]  ^2  S^n  k2X i , 

[^0406  ]  k-2  C0S  k-2X\ 

[^o5oi  ]  coskjXj , 

[^0502  ]  sinkj.Xj, 

[^0505  ]  k2  cos  k2xx , 

[^()6o i  ]  sin^x,,  [ 

^0602]  ^1  coskjXj , 

[^0605  ]  ^-2  sink2Xj 

[^0705  ]= cos  Ar2xc2  > 

[-^0706  ]=  s^n  k2x2 , 

[^0707  ]  =  X2  5 

[^0709  ]  ~  COS  kyX2  , 

[^0710  ] =  sin  k2x2 ,  | 

.^071 1  ]  =  X2  5 

[-^0805 \  —  —  k2S:mk2X2, 

[^0806  ]  ^2  COSk2X2  ,  | 

[^0807  ]=1  ’ 

|7o81o]  —  _  ^3  COS  k2X2  , 

1 

i  II 

OO 

[^0905  ]  ^"2  COSk2X. 

[^0909 ]  ^3  COSk3X2  , 

[^0910 ]  k^  sin  k2x2 , 

[^1005  ]  ^2  sink2x2 , 

[^1009  ]  ^3  sin  k2x2 , 

faJ-*3  COS^3X2  ’ 

|  jo9  ]  cosk3x3, 

[^1 1 1 1  ]  =  *^3  ’ 

fcmH 

[7)209  ]  =  —  ^3  C0S  ^3X 

fcno]  sin  A^3x3 ,  [^1 111]  x3’  [-^1112]  [^1209 ]  £3  cos k3x3, 

[^1210  ]~  ^3  sinA:3x3 

The  others  of  [7^2ri2]  are  all  zero. 

For  a  nontrivial  solution  of  Equation  (18),  the  determinant  of  the  above  matrix  is  set  equal 
to  zero,  \Tl2xl2\  =  0. 

To  find  a  particular  solution,  we  take 

y0  =  Y0eirx,  y  =  Be'rx ,  i  =  imaginary  number.  (19) 

Substitution  Equation  (19)  into  Equation  (3)  gives 


1  3 

=  —  ^  Y0J  (cosh  kjX  -  sinh  kfx) 
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Therefore  the  general  solution  is 


3  i 

yg=Y[cn  cos  kjX+  Cj2  sin  kjX  +  Cj3x  +  C/4  —  70 ;  (cosh  ky  -  sinh  kfx)]  (2 1 ) 


7=1 


The  equations  of  motion  for  dynamic  buckling  is  the  same  Equation  (1)  that  was 
mentioned  above.  After  dividing  through  by  El,  we  assume  parameters  as  follows. 


2  P  i  I  2  P 

k  =  — ,  r  =  — ,  c  =  — 
El  A  p 


where  A  is  an  area  and  p  is  a  mass  density  of  beam. 
Then  the  Equation  ( 1 )  becomes 


(22) 


3  y  ,  12  d  y  1  d  y  _  ,^2  d  y0 


dx4 


-  +  kz^y  +  - 


8xz  r2c 2  dt 2 


=  -kz 


o>x2 


(23) 


In  case  of  a  stepped  beam,  the  equation  of  motion  is 

54v;.  ,  2  d2yj  .  1  d2y_j  2  d2y0j 


J  \  kA  J  \  J  ~  k  _ 

4  A  o  2  ^  „  2  2  r\.2  Kj  ^2 


dx 4  2  dx2  r  c  dt  J  dx 

J  J 

We  assume  the  solution  of  Equation  (24)  as  follows. 
y(x,t)  =  f(x)g(t ) ,  v0(x,0)  =  f0(x) 

where, 

3 


j=l,2,3 


(24) 


(25) 


3  1 

f{x)  =  y{x)  =  YJ[Cn  cos  kjX  +  Cj2  sin  kjX  +  Cj3x  +  Cj4  -  —  Y0J  (cosh  kjX  -  sinh  kjX )]  (26) 


7=1 


3  3 

-k,: 


fo  O)  =  To  (x)  =  Z  Y0 je'V  =  Z  Y0je  JX  =  Z  Y0j  (cosh  krx  -  sinh  kj x)  (27) 

7=1  7=1  7=1 

Substituting  Equation  (25)  into  Equation  (24)  and  arranging,  the  following  equation  is 
given. 


~H2HA)g(t)-HlH2HAg(t)  =  -HXH2H4 
Dots  indicate  differentiation  with  respect  to  time.  Where, 


(28) 


Hx  =  k) 


(29) 
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(30) 


H2  =  Y0j 


fjC. 


H4  =  cosh  kjX  -  sinh  kjX 


II 5  =  C cos  k^  +  CJ2  sin  k  .x  +  Cj  Vx  +  CJ4 


Rewriting  Equation  (28) 

g(t)  +  Sg(t)  =  S 

where, 

s  =  -^  =  - 

5. 


-hxh2h4 


H2(Hs-^H2H4) 


(31) 

(32) 

(33) 

(34) 

(35) 


The  homogeneous  solution  of  dynamic  buckling  can  be  expressed  as  follows.  Substituting 
g(t)=ei'  into  Equation  (34),  a  characteristic  equation  is  A2  =  -S  . 

In  case  of  S  >0;  A  =  ±y[si  =  ±Si 


gh  ( t )  =  Wj  cos Sst  +  u2  sin  Sst  (36) 

In  case  of  S  <  0;T  =  ±Sr 

gh  ( t )  =  uxeSr‘  +  u2e~s,t  =  Mj  cosh  Srt  +  ux  sinh  Srt  +  u2  cosh  Srt  -  u2  sinh  Srt 


=  (u4  +w2)coshS',.t  +  (iq  -u2) sinh Srt  (37) 

In  case  of  5*  =  0;A  =  0 

gh(t)  =  ui+u2t  (38) 

The  particular  solution  of  dynamic  buckling  is  as  follows.  Substituting  g(t )  =  A  into 
Equation  (34),  a  characteristic  equation  is  A=  1. 

gP(t)  =  1  (39) 

Therefore  the  general  solution  is 

gg(t)  =  gh(t)  +  gP(t)  (4°) 
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y(x,t)  =  f(x)gg(t)  =  f(x)g(t)  (41) 

The  beam  is  assumed  to  be  initially  at  rest.  Also,  y  is  measured  from  the  initial 
displacement  Y0 ,  so  the  initial  conditions  are 


v(x,0)  =  y(x,0)  =  0  (42) 

Applying  these  to  Equation  (41)  yields  u]  and  u2  for  each  case.  The  final  solutions  are  then 

3  i 

y(x,t)  =  2K,  cos  kjX  +  Cj2  sin  kjX  +  CJ3x  +  Cj4  —  Y0j  (cosh  kx  -  sinh  kjX)\ 


j= i 


x(-cos5'st  +  l) ,  ^>0 


(43) 


3 

y(x,t)  =  sx,  cos  kjX  +  CJ2  sin  kjX  +  Cj3x  +  CjA 
j= i 

x(-cosh5',f +  1) ,  ,S<0 
y(x,t)  =  0  ,  S  =  0 


1 

2 


Y0j  (cosh  kjX  -  sinh  Ayx)] 


(44) 

(45) 


B.  A  CANTILEVER  BEAM 

If  the  external  force  is  applied  at  the  free  edge  of  a  cantilever  beam  with  stiffness 
variation  as  in  Figure  2,  the  equations  of  motion  are  given  as 


t 

H 


Figure  2.  Stepped  Cantilever  Beam  Subjected  to  a  Vertical  Load 
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d2v  . 

J7  J  '  '  J  i  ~  A  _ LJ— 


d4yj 


+pA 


=  0,  j=l,2,3 


'rj  dx4  1  1  dt2 

We  assume  that  the  solution  of  Equation  (46)  is 

yj=Yj(x)eM 

Substituting  Equation.  (47)  into  Equation  (46) 

8%  2 

YJ{x)  =  YJ=Cjei* 

Substituting  Equation  (49)  into  Equation  (48) 

(EjIjA4  -  PjAjCo2  )Cj=0 
For  the  nontrivial  solution  of  Equation  (50) 

\Ejl/r  PjAjCo1^ 

We  can  get  twelve  roots  of  A  from  Equation  (51).  Finally,  we  have 
W'ECj/''  .  j=l .2,3 

n= 1 

The  boundary  and  the  continuity  condition  to  get  C  .  of  Equation  (52)  are 
At  x  =  0 

Ti(0)  =  0 

y[(0)  =  o 
At  x  =  x,,x2 

yl(xl)  =  y2(xl) 
y[(xl)=y'2(xl) 
y"(x])  =  y2(xt) 
v;"(x1)=v;"(x1) 

y  2  ("^2 ) = T3  (x2 ) 

T2  (x2)= y'3(x2) 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
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y"(x2)  =  y”(x2) 
y'2(x2)=y^x2) 


At  x  =  x. 


y''(-h)  =  () 

y2(x3)  =  -  F° 


From  Equation  (53  through  64),  we  can  get  12  C .  as  follows. 


r  = 

^34 


Fn 


jA, /,  /?,  cos (f?3X3) 


r 


sin(f?3x3) 

cos(f?3x3) 


C 


34 


sin(7f.x2 )  cos(i?3x2) 

cosh(f?3x2)  33  cosh(f?3x2)  34 


(61) 

(62) 

(63) 

(64) 


sin(i?3x2)  c  |  cos(i?3x2)  c  |  sin (R3x2)  c 
cosh(f?3x2)  32  cosh(f?3x2)  33  cosh(f?3x2) 


C 


24 


c 


23 


R3  sinh(7?3x2)  n 

'-'31 

R2  cos(f?2x2) 
sin(i?2x2)  c  + 
cos(f?2x2)  24 


+  R3  cosh(/?3x2)  f?3  sin(/?3x2 )  ^  +  R3  cos (R3x2) 

R2cos(R2x2)  32  R2cos(R2x2)  33  R2cos(R2x2)  34 
cosh(/?3x2)  ^  +  sinh(/?3x2)  ^  +  cos (R3x2)  ^  +  sin(f?3x2) 
cos(f?2x2)  31  cos(f?2x2)  32  cos(f?2x2)  33  cos(R2x2) 


C 


34 


si n(i?2Xi)  c  |  cos(i?2Xj)  c 
cosh(f?2Xj)  23  cosh(f?2Xj) 


sinh^x,)  |  cos(i?2X!)  |  sin(/?2x, ) 
cosh(f?2Xj)  22  cosh(f?2Xj)  23  cosh(f?2Xj)  24 

f?2  sinh(f?2Xj)  +  R2  cosh(f?2Xj)  f?2  sin(f?2Xj)  +  f?2  cos(f?2x1) 
f?j  cos(f?1x1)  21  f?[  cos(f?jXj)  22  f?,  cos(f?1x1)  23  f?,  COS(f?,Xj) 
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sin  (Rft)  cosh  {R2xx)  smli^x,)  cos  (R2x{) 

13  COS(i?jXj)  14  COS(i?jXj)  21  COS(i?jXj )  22  COS(i?,Xj) 

C12  =  -c14 

cn=-c13 


|  sin(/?2x, )  c 

COS^jXj) 
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II.  NUMERICAL  ANALYSIS  AND  CONSIDERATIONS 


For  investigating  what  the  stiffness  variation  of  structures  affect  on  a  system  under 
an  impact  load,  a  numerical  analysis  is  carried  out  for  several  models.  The  MSC/PATRAN 
[Ref.  11]  for  finite  element  modeling  and  the  MSC/NASTRAN  program  [Ref.  12]  for 
transient  analysis  are  used.  In  this  report,  a  cantilever  beam  subjected  to  a  vertical  load  as 
shown  in  Figure  2  and  a  simply  supported  beam  subjected  to  an  axial  load  as  in  Figure  3 
are  analyzed.  Then  a  transient  analysis  is  carried  out  also  for  a  cylindrical  shell  of  a 
cantilever  condition.  The  purpose  is  to  confirm  whether  the  analysis  result  of  a  beam  of  one 
dimensional  configuration  coincides  with  that  of  a  cylindrical  shell  of  three  dimensional 
configuration  such  as  Figure  4  or  not.  In  a  preliminary  design  stage,  to  model  the  large 
structures  of  a  cylindrical  shell  type  like  ship,  missile  and  airplane  as  one  dimensional  beam 
is  very  important  to  save  time  and  effort.  The  response  of  structures  with  the  frequency 
variation  of  impact  loads  are  also  investigated  for  the  beam  and  the  cylinder  of  a  cantilever 
condition.  Table  1  shows  the  load  and  the  boundary  conditions  applied  in  transient 
analyses.  Table  2  shows  the  dimension  and  the  physical  property  of  material  of  each  model 
used  for  a  numerical  analysis. 
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Figure  3.  Stepped  Simply  Supported  Beam  Subjected  to  an  Axial  Load 
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Figure  4. 


Configuration  of  a  Cylindrical  Shell 


Table  1.  Applied  Loads  and  Boundary  Conditions  for  Transient  Analysis 


Load 

and 

Type 

Load:  F(t)  =  F0  sin  cot 

Boundary 
Conditions 
and  Models 

Shape 

Direction 

Fo(kgf) 

co(  hz) 

Time 
(t ,  sec) 

ti 

t2 

Load 

1 

Half  Sine 
Wave 

Axial 

161,167 

50 

0.0 

0.0100 

Simply 

Supported, 

Beam(ssb) 

Load 

2 

Half  Sine 
Wave 

Vertical 

1,600 

50 

0.0 

0.0100 

Cantilever, 
Beam  and 
Cylindrical 
Shell(cb/ccs) 

Load 

3 

Half  Sine 
Wave 

Vertical 

1,600 

100 

0.0 

0.0050 

Cantilever, 
Beam  and 
Cylindrical 
Shell(cb/ccs) 

Load 

4 

Half  Sine 
Wave 

Vertical 

1,600 

200 

0.0 

0.0025 

Cantilever, 
Beam  and 
Cylindrical 
Shell(cb/ccs) 
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Table  2.  Dimension  and  Physical  Data  of  Beam  and  Cylinder  Model 

(Unit:  mm) 


Dimension 

and 

Model  Type 

Section  1 

Section  2 

Section  3 

O.D* 

I.D*  Fength 

O.D* 

I.D* 

Length 

O.D* 

I.D* 

Length 

Model  1 

302 

298  1000 

302 

298 

1000 

302 

298 

1000 

Model  2 

302 

298  1000 

302 

299 

1000 

302 

298 

1000 

Model  3 

302 

298  1000 

302 

300 

1000 

302 

298 

1000 

Model  4 

302 

298  1000 

302 

301 

1000 

302 

298 

1000 

Material  A1 2024-t3 

Elastic  Modulus  7,380  (kgf/mnr) 

Mass  Density  2.853E-10  (kgf.sec2/mm4) 

Poisson’s  Ratio  0.33 

*  :  O.D  and  I.D  are  acronym  for  Outer  and  Inner  Diameter. 

A.  A  SIMPLY  SUPPORTED  BEAM  UNDER  AN  AXIAL  IMPACT  LOAD 

1.  The  Effect  of  a  Stiffness  Variation 

Figure  5  through  14  show  the  displacements  and  the  stresses  of  a  simply  supported 
beam  subjected  to  an  axial  impact  load,  the  load  1  in  Table  1.  The  purpose  is  to  understand 
the  phenomenon  varying  the  displacements  and  the  stresses  by  the  thickness  reduction  of 
section  2  in  Table  2,  namely,  the  stiffness  change  of  structures. 

Figure  5,  which  is  in  case  of  model  1  without  the  thickness  variation  of  each 
section,  shows  axial  displacements  on  a  time  domain.  From  the  graphs,  we  know  that  the 
axial  displacements  are  the  largest  at  the  sta.3000  applied  an  axial  impact  and  the  others  are 
decreasing  in  order  of  sta.  2000,  1500,  1000  and  500.  The  displacements  have  the 
characteristics  appearing  mainly  for  0.01  sec  applying  an  axial  impact,  then  reducing 
abruptly  and  disappearing  finally. 
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Figure  6  is  the  stress  curve  of  model  1.  The  stresses  are  the  largest  at  the  sta.  1500, 
the  mid-position  of  the  section  2  among  three  sections.  It  is  natural  that  the  largest 
compressive  stress  of  a  simply  supported  beam  subjected  to  an  axial  impact  is  to  come  out 
at  the  mid  position  of  the  beam.  The  stress  like  the  displacement  happened  mainly  for 
applying  an  axial  impact,  and  it  disappeared.  The  magnitude  of  the  axial  impact  is  161,167 
kg,  a  critical  buckling  load  of  the  model  1.  Therefore,  from  the  figure,  we  know  that  the 
stresses  of  all  models  exceed  the  ultimate  tensile  strength  of  A12024-t3,  45  kg/mm  , 
already. 

Figure  7  is  the  displacement  curve  of  model  2  that  the  thickness  of  section  2  is 
1.5mm.  A  general  trend  is  similar  to  Figure  5,  but  the  magnitudes  of  the  displacements  are 
increasing  large. 


16 


Figure  6.  Stresses  of  Model  l(ssb,  2mm)  under  Load  1  (axial,  50hz) 


Figure  7.  Displacements  of  Model  2(ssb,  1.5mm)  under  Load  1  (axial,  50hz) 
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Figure  8  is  the  stress  curve  of  model  2.  The  compressive  stress  is  the  maximum  at 
the  mid  position  of  the  beam  and  increased  than  the  Figure  6.  It  is  because  of  a  stiffness 
weakening  by  the  thickness  reduction  of  section  2. 

Figure  9  and  10  are  the  displacement  and  the  stress  graphs  of  model  3.  The 
thickness  of  section  2  is  1 .0mm  which  is  the  half  of  model  1 .  The  trend  of  the  displacement 
and  the  stress  is  similar  to  Figure  7  and  8,  but  the  magnitude  increases. 

Figure  1 1  and  12  are  of  model  4  that  the  thickness  of  section  2  is  0.5mm.  It  shows 
that  the  displacement  is  about  double  and  the  stress  is  about  4  times  larger  than  Figure  5 
and  6  of  model  1 . 

From  Figure  5  through  12,  we  know  that  the  displacement  and  the  stress,  at  the  sta. 
3000  and  at  the  sta.  1500,  are  increased  due  to  the  stiffness  weakening  of  structures, 
namely,  the  thickness  reduction  of  section  2. 


Figure  8.  Stresses  of  Model  2(ssb,  1.5mm)  under  Load  l(axial,  50hz) 
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Stress(Kg/mm  )  £  X  displacement(mm) 


Displacements  of  model  3(ssb,  1mm)  under  load  l(axial,  50hz) 


Figure  10.  Stresses  of  Model  3(ssb,  hum)  under  Load  l(axial,  50hz) 
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Figure  11.  Displacements  of  Model  4(ssb,  0.5mm)  under  Load  1  (axial,  50hz) 


Figure  12.  Stresses  of  Model  4(ssb,  0.5mm)  under  Load  l(axial,  50hz) 
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Figure  13  is  the  graph  that  done  a  comparative  study  on  the  displacements  of 
models  at  the  sta.  3000  occurred  the  largest  displacement.  The  increasing  rate  of  the 
displacement  is  larger  than  the  decreasing  rate  of  stiffness,  and  the  order  is  model  4,  3,  2 
and  1 .  The  time  occurring  the  largest  displacement  is  delaying  by  the  thickness  reduction  of 
the  section  2. 

Figure  14  is  to  compare  the  stresses  of  models  at  the  sta.  1500  occurred  the  largest 
stress.  The  stress  has  also  a  similar  trend  like  the  displacement,  and  the  stress  of  model  4  is 
the  largest. 

Table  3  is  the  results  gotten  by  a  normal  mode  analysis  and  a  transient  analysis  to 
investigate  what  the  thickness  reduction  of  section  2  affects  on  the  dynamic  characteristics 
of  structures.  The  area  moments  of  inertia  and  the  first  natural  frequency  are  reduced  by  the 
thickness  reduction  of  models,  but  the  displacements  and  the  stresses  are  increased  about 
double  and  four  times  each. 


Figure  13.  Displacements  of  Models(ssb)  at  sta.  3000  under  Load  1  (axial,  50hz) 
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Figure  14.  Stress  of  Models(ssb)  at  sta.  1500  under  Load  l(axial,  50hz) 


Table  3.  Analysis  Results  of  Simply  Supported  Beams 


Results 

and 

Model 

Type 

Area  Moment 
of  Inertia 

First  Natural 
Frequencies 

Transient  Analysis 

Mag.*1 

(mm4) 

Rate 

(%) 

Mag.*1 

(hz) 

Rate 

(%) 

Maximum 

Displ.*2 

(mm) 

Rate 

(%) 

Maximum 

Stress 

2 

(kgf/mm  ) 

Rate 

(%) 

Model  1 

2.12E07 

100 

91.4 

100 

33.0 

100 

107 

100 

Model  2 

1.60E07 

75 

91.0 

99.6 

37.2 

113 

139 

130 

Model  3 

1.07E07 

50 

88.2 

96.5 

46.1 

140 

205 

192 

Model  4 

5.38E06 

25 

77.6 

84.9 

73.8 

224 

405 

379 

Mag.  is  an  abbreviation  of  magnitude. 

2 

*  :  Displ.  is  an  abbreviation  of  displacement. 
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B.  A  CANTILEVER  BEAM  UNDER  A  VERTICAL  IMPACT  LOAD 

1.  The  Effect  of  a  Stiffness  Variation 

From  Figure  5  through  14,  we  knew  that  the  stiffness  variation  of  a  simply 
supported  beam  subjected  to  axial  impact  load  affects  large  on  the  displacement  and  the 
stress  of  structures.  Therefore,  we  are  going  to  investigate  that  how  a  vertical  impact  affects 
on  a  cantilever  beam,  one  among  the  boundary  conditions  used  often  in  a  structural  design. 
A  vertical  load  is  selected  because  the  bending  stress  often  induces  the  maximum  stress  in 
structure.  The  magnitude  of  an  impact  load,  1600kg,  was  determined  by  a  static  analysis. 
That  is  the  load  just  before  the  working  stress  is  over  the  allowable  stress  of  material. 

Figure  15  is  the  displacement  graph  of  model  1.  The  displacement  is  increasing  in 
proportion  to  the  arm  length  of  a  bending  moment  because  the  stiffness  of  model  1  is  not 
changed.  The  difference  from  a  simply  supported  beam  subjected  to  an  axial  impact  is  that 
the  displacement  is  decreasing  a  little  and  occurring  alternately  after  0.01  sec  done  an 
impact.  It  means  that  a  cantilever  beam  vibrates  above  and  below  because  the  vertical 
impact  load  is  applied  at  the  free  edge  of  that. 

Figure  16  is  the  stress  curve  of  model  1.  Also,  the  stress  is  occurring  at  all  position 
of  a  cantilever  beam  after  0.01  sec  done  an  impact.  This  is  also  a  different  point  from  a 
simply  supported  beam  under  axial  load.  The  magnitude  order  of  stress  is  opposite  from  the 
displacement.  It  is  because  a  bending  stress  is  the  largest  at  the  longest  arm  of  a  bending 
moment  as  the  static  analysis. 

Figure  17  is  the  displacement  curve  of  model  2  that  the  thickness  of  section  2  is 
1.5mm.  The  displacement  is  increasing  than  model  1,  but  the  trend  is  similar  to  model  1. 

Figure  18  is  the  stress  curve  of  model  2.  The  stress  is  increasing  and  the  trend  is 
similar  to  model  1 . 

Figure  19  is  to  show  the  displacements  of  model  3  that  the  thickness  of  section  2  is 
1.0mm.  The  general  trend  is  similar  to  model  2,  but  the  magnitudes  of  displacements  are 
increasing. 
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Figure  15.  Displacements  of  Model  l(cb,  2mm)  under  Load  2(vertical,  50hz) 


Figure  16.  Stresses  of  Model  l(cb,  2mm)  under  Load  2(vertical,  50hz) 
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Figure  17.  Displacements  of  Model  2(cb,  1.5mm)  under  Load  2(vertical,  50hz) 


Figure  18.  Stresses  of  Model  2(cb,  1.5mm)  under  Load  2(vertical,  50hz) 
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Figure  19.  Displacements  of  Model  3(cb,  1mm)  under  Load  2(vertical,  50hz) 


Figure  20  shows  the  stress  of  model  3.  The  different  from  the  stress  of  model  2  is 
that  the  stress  is  increasing  abruptly  at  the  sta.  1000  and  1500.  The  stress  distribution  is 
changing  by  that  the  stiffness  of  section  2  is  reduced  to  50%. 

Figure  21  is  the  displacement  curve  of  model  4  that  the  thickness  of  section  2  is 
0.5mm.  The  displacement  is  the  largest  among  models.  But  the  trend  is  also  similar  to 
others. 

Figure  22  is  the  stress  curve  of  model  4.  It  shows  that  the  stress  is  the  largest  at  the 
sta.  1500  of  the  mid  position  of  a  cantilever  beam.  It  is  because  the  stiffness  reduced  to 
75%  than  that  of  model  1,  and  the  stress  values  are  exceeding  the  ultimate  tensile  strength 
of  an  A12024-t3,  45kg/mm2. 

Table  4  is  the  analysis  results  of  cantilever  beams  subjected  to  a  vertical  impact  of 
load  2.  The  stress  and  the  displacement  are  increased  by  thickness  reduction  of  section  2. 
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Figure  20.  Stresses  of  Model  3(cb,  1mm)  under  Load  2(vertical,  50hz) 


Figure  2 1 .  Displacement  of  Model  4(cb,  0.5mm)  under  Load  2(vertical,  50hz) 
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Figure  22.  Stresses  of  Model  4(cb,  0.5mm)  under  Load  2(vertical,  50hz) 


Table  4.  Analysis  Results  of  Cantilever  Beams 


Model 

Type 

and 

Results 

Area  Moment 
of  Inertia 

First  Natural 
Frequencies 

Transient  Analysis 

Mag.*1 

(mm4) 

Rate 

(%) 

Mag.*1 

(hz) 

Rate 

(%) 

Maximum 

Disp.*2 

(mm) 

Rate 

(%) 

Maximum 

stress 
(kgf/mrn  ) 

Rate 

(%) 

Model  1 

2.12E07 

100 

33.0 

100 

101.8 

100 

33.8 

100 

Model  2 

1.60E07 

75 

32.8 

99.4 

110.5 

108. 

5 

34.0 

100.6 

Model  3 

1.07E07 

50 

31.7 

96.1 

125.4 

123. 

2 

36.3 

107.4 

Model  4 

5.38E06 

25 

28.0 

84.8 

160.9 

158. 

1 

62.3 

184.3 

*2 


:  Mag.  is  an  abbreviation  of  magnitude. 

:  Displ.  is  an  abbreviation  of  displacement. 
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It  is  investigated  the  displacements  and  the  stresses  in  all  cantilever  models  with  a 
stiffness  variation  subjected  to  a  vertical  impact  load.  The  displacement  is  the  largest  at  the 
free  edge  of  a  cantilever  beam  and  the  farther  from  that,  it  is  lessened.  In  case  of  stress,  it  is 
the  largest  at  the  fixed  position  of  a  cantilever  beam  but  it  shows  the  phenomenon  that  the 
stress  is  moving  to  the  place  that  the  stiffness  reduction  is  large. 

2.  The  Effect  of  a  Frequency  Variation 

Figure  23  through  28  are  to  investigate  what  the  frequency  change  of  an  impact  load 
affects  on  the  models 

Figure  23  shows  the  displacements  of  models  at  the  sta.  3000  when  the  impact  of  50 
hz  applies  a  cantilever  beam.  The  displacement  of  model  4  is  the  largest,  and  the  response 
is  delaying  than  model  1,  2  and  3.  The  reason  is  that  the  stiffness  reduction  rate  of  section  2 
is  larger  than  other  models. 

Figure  24  is  the  stress  curve  at  the  sta.  1500  which  is  the  position  occurring  the 
maximum  stress  when  the  impact  of  50  hz  is  applied.  The  stress  of  model  4  is  larger  than 
those  of  model  1,  2  and  3.  This  is  one  of  important  factors  which  should  be  considered 
when  design  structures.  Model  1,  2  and  3  are  safe  because  the  stresses  of  those  are  smaller 
than  the  ultimate  tensile  strength  of  material,  but  model  4  is  not  safe  for  it  is  over  he 
ultimate  tensile  strength. 

Figure  25  is  the  displacement  curves  of  models  subjected  to  the  impact  load  of 
lOOhz.  The  displacements  are  decreasing  than  the  models  of  50hz.  The  reason  is  that  the 
applied  energy  is  reduced  and  it  makes  shorten  the  response  of  the  structures  because  the 
magnitude  of  an  impact  does  not  change  but  time  only  shortens. 

Figure  26  is  the  stress  curves  of  models  subjected  to  an  impact  load  of  lOOhz.  The 
stresses  like  the  displacement  are  decreasing  than  those  of  the  50hz  models.  And  another 
difference  from  the  results  of  50hz  is  that  the  stress  value  is  the  maximum  at  the  second 
peak  than  at  the  first  one.  Because  it  is  added  before  the  first  peak  does  not  finish 
completely. 
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Figure  23.  Displacements  of  Models(cb)  at  sta.  3000  under  Load  2(vertical,  50hz) 


Figure  24.  Stresses  of  Models(cb)  at  sta.  1500  under  Load  2(vertical,  50hz) 
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Figure  25.  Displacements  of  Models(cb)  at  sta.  3000  under  Load  3(vertical,  lOOhz) 


Figure  26.  Stresses  of  Models(cb)  at  sta.  1500  under  Load  3 (vertical,  lOOhz) 
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Figure  27  is  the  displacement  curves  of  models  under  an  impact  load  of  200hz.  The 
displacements  are  decreased  large  because  the  shock  energy  is  reduced  by  the  reduction  of 
an  applying  time. 

Figure  28  is  the  stress  curves  of  models  under  an  impact  load  of  200hz.  The  stresses 
are  decreasing  than  the  models  of  lOOhz,  but  the  maximum  stress  value  is  the  second  peak 
as  that  of  lOOhz. 

It  has  been  investigated  what  the  frequency  variation  of  impact  load  affects  on  the 
stresses  of  models.  The  increment  of  a  frequency  only  of  an  applied  impact  load  would 
rather  reduce  the  stress  of  each  model  because  of  the  reduction  of  an  applied  energy. 


Figure  27.  Displacement  of  Models(cb)  at  sta.  3000  under  Load  4(vertical,  200hz) 
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Figure  28.  Stresses  of  Models(cb)  at  sta.  1500  under  Load  4(vertical,  200hz) 


C.  A  CYLINDRICAL  SHELL  UNDER  A  VERTICAL  IMPACT  LOAD 

Till  now,  studies  were  conducted  on  the  simply  supported  beam  and  the  cantilever 
beam  subjected  to  an  axial  and  a  vertical  impact  load.  We  are  going  to  apply  the  same  load 
that  applied  a  cantilever  beam  to  a  cylindrical  shell  as  in  Figure  4,  and  to  compare  the 
results  of  a  beam  and  a  cylindrical  shell.  Physical  properties  of  material,  force  and 
boundary  condition  used  on  this  occasion  are  the  same.  In  the  preliminary  design  stage, 
modeling  and  analyzing  three  dimensional  cylinder  to  one  dimensional  beam  are  important 
in  the  point  of  view  saving  time  and  effort. 

Figure  29  and  30  are  the  displacement  and  the  stress  curves  of  the  cylindrical  shell 
with  a  cantilever  boundary  condition  subjected  to  a  vertical  impact 

From  Figure  29,  the  displacements  of  a  cylindrical  shell  are  the  maximum  at  the  sta. 
3000,  and  the  order  is  sta.  2500,  2000,  1500  and  500.  It  shows  the  same  trends  as  a 
cantilever  beam. 
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Figure  29. 


Displacements  of  Model  l(ccs,  2mm)  under  Load  2(vertical,  50hz) 


Figure  30.  Stresses  of  Model  l(ccs,  2mm)  under  Load  2(vertical,  50hz) 
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In  case  of  stress,  the  order  of  the  magnitude  is  sta.  500,  1000,  1500,  2000,  2500  and 
3000.  This  coincides  the  magnitude  order  of  a  bending  moment  and  the  trend  is  the  same  as 
beam. 

Figure  31  through  42  is  the  graphs  to  compare  the  analysis  results  of  cantilever 
beams  and  cylindrical  shells  with  frequency  variation  at  the  position  of  maximum  stress. 

Figure  31  is  the  stress  curves  of  model  1  on  the  frequency  50hz  at  the  sta.  1500.  It 
shows  a  similar  trend,  but  the  stress  of  a  cylindrical  shell  is  larger  than  that  of  beam.  It  is 
because  of  the  model  characteristics  of  a  beam  and  a  cylindrical  shell.  It  is  the  result  to 
compare  the  mean  stress  of  the  element  for  a  beam  with  that  at  the  farthest  of  the  element 
for  a  cylindrical  shell.  The  stiffness  of  cylindrical  shell  is  reduced  than  that  of  beam 
because  the  finite  element  of  a  cylindrical  shell  is  not  fine  mesh.  Whether  the  element  of  a 
cylindrical  shell  make  a  fine  mesh  or  not,  that  is  the  best  way  to  decide  when  required  more 
detail  analysis  after  an  analysis  of  a  beam  model.  This  is  to  reduce  a  computing  time  and  an 
effort.  The  time  occurring  the  peak  value  of  a  cylindrical  shell  is  delaying  than  that  of  a 
beam  after  the  third  peak  value.  It  takes  time  to  respond  a  structure  after  an  external  load  is 
applied  because  a  cylindrical  shell  is  three  dimensional  structure  and  is  more  complex  than 
a  beam. 

Figure  32,  33  and  34  are  the  result  applied  the  load  of  50hz  for  the  model  2,  3  and  4. 
It  shows  a  similar  trend  with  model  1.  But  the  stresses  are  increasing  little  by  little.  It  is 
because  of  the  stiffness  reduction  of  section  2. 

Figure  35,  36,  37  and  38  are  the  stress  curves  of  the  case  applied  the  load  of  100  hz 
for  each  model.  The  analysis  results  for  a  beam  and  a  cylindrical  shell  are  similar,  and  the 
less  the  stiffness,  the  more  the  stress.  But  the  stresses  are  decreasing  than  that  of  50hz.  It  is 
because  the  application  time  is  shortened  and  the  energy  is  lessened.  And  the  second  peak 
is  the  largest  because  it  is  added  before  the  first  peak  does  not  finish. 

Figure  39  through  42  are  the  case  applied  the  load  of  200hz  for  each  model.  It 
shows  a  similar  trend  with  that  of  lOOhz.  But  the  stresses  are  decreasing  than  that  of  lOOhz. 
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Stress(Kg/mm  )  r"  Stress(Kg/mm  ) 


Stresses  of  Model  l(cb/ccs,  2mm)  at  sta.  1500  under  Load  2(vertical,  50hz) 


Figure  32.  Stresses  of  Model  2(cb/ccs,  1.5mm)  at  sta.  1500  under  Load  2(vertical, 

50hz) 
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Stress(Kg/mm  )  Stress(Kg/mm  ) 


Stresses  of  Model  3(cb/ccs,  1mm)  at  sta.  1500  under  Load  2(vertical,  50hz) 


Figure  34.  Stresses  of  Model  4(cb/ccs,  0.5mm)  at  sta.  1500  under  Load  2(vertical, 

50hz) 


37 


Figure  35.  Stresses  of  Model  l(cb/ccs,  2mm)  at  sta.  1500  under  Load  3(vertical,  lOOhz) 


Figure  36.  Stresses  of  Model  2(cb/ccs,  1.5mm)  at  sta.  1500  under  Load  3(vertical, 

lOOhz) 
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Stress(Kg/mm  )  •  Stress(Kg/mm  ) 


Stresses  of  Model  3(cb/ccs,  1mm)  at  sta.  1500  under  Load  3(vertical,  lOOhz) 


Figure  38.  Stresses  of  Model  4(cb/ccs,  0.5mm)  at  sta.  1500  under  Load  3(vertical, 

lOOhz) 


39 


Stress(Kg/mm  )  1°  Stress(Kg/mm  ) 


Stresses  of  model  l(cb/ccs,  2mm)  at  sta.  1500  under  load  4(vertical,  200hz) 


Figure  40.  Stresses  of  model  2(cb/ccs,  1.5mm)  at  sta.  1500  under  load  4(vertical, 

200hz) 
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Figure  41.  Stresses  of  model  3(cb/ccs,  1mm)  at  sta.  1500  under  load  4(vertical,  200hz) 


Figure  42.  Stresses  of  model  4(cb/ccs,  0.5mm)  at  sta.  1500  under  load  4(vertical, 

200hz) 
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It  has  been  investigated  on  the  difference  between  a  beam  model  and  a  cylindrical 
shell  model.  The  results  were  similar  if  the  condition  was  the  same. 
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IV.  CONCLUSIONS 


The  research  results  on  the  dynamic  response  of  stepped  beams  and  stepped 
cylindrical  shells  subjected  to  an  impact  load  are  as  follows. 

1 .  The  theoretical  studies  on  the  analytic  solutions  of  a  simply  supported  beam  and  a 
cantilever  beam  were  carried  out. 

2.  It  was  analyzed  that  a  stiffness  reduction  affects  on  the  dynamic  characteristics  of 
a  simply  supported  beam  subjected  to  an  axial  impact.  By  the  thickness  reduction  of 
models,  the  displacement  and  the  stress  were  increased  about  double  and  four  times  each. 

3.  In  case  of  a  cantilever  beam  subjected  to  a  vertical  impact  load,  the  structural 
response  was  increased  by  the  stiffness  reduction  as  the  simply  supported  beam.  The  stress 
was  the  largest  at  the  fixed  position  of  beam,  but  it  moved  and  increased  toward  which  the 
stiffness  reduction  was  large. 

4.  The  analysis  results  of  a  beam  and  a  cylindrical  shell  were  similar  if  the  physical 
property  of  a  material,  boundary  and  load  conditions  were  same. 

5.  The  stiffness  reduction  was  a  factor  increasing  a  stress  and  a  displacement  of 
structures  regardless  of  frequency  variation.  Because  the  magnitude  of  an  impact  load  does 
not  change  but  time  only  shortens. 

It  is  still  needed  parametric  study  as  an  initial  deflection  and  an  eccentric  load.  But, 
we  hope  this  is  to  be  a  little  help  and  is  to  be  the  first  step  to  reduce  an  effort  and  time  in  the 
preliminary  design  stage  of  complex  structures  like  ship,  airplane  and  missile. 
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